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Abstract 

Nucleation of vortices in rotating superfluid by spin-up and rapid ther- 
mal quench is discussed in the framework of the time-dependent Ginzburg- 
Landau equation (TDGLE). An analysis of the instability in inhomogeneous 
rotationally-invariant system results in the expression for the critical rota- 
tional velocity. A stability analysis of multicharged vortices is presented. It 
is shown that they are very long-living objects with lifetime inversely pro- 
portional to the dissipation rate. It was found by numerical and analytical 
solution of the TDGLE that vortex nucleation by rapid thermal quench in the 
presence of superflow is dominated by a transverse instability of the moving 
normal-super fluid interface. 
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I. INTRODUCTION. 



There are two major ways to generate vorticity in rotating superfluid: to spin-up a bucket 
with superfluid helium from initially steady state without any vortices to a state with a ro- 
tation velocity Q > Q c , where Q c is the critical rotation speed for a vortex nucleation, and to 
quench it thermally from temperature, T, above the superfluid transition temperature, T\, 
to temperature below it. These two different scenaria can be described by the same model, 
namely time-dependent Ginzburg-Landau equation (TDGLE) for the superfluid order pa- 
rameter corrected for presence of normal component. This equation for superfluid helium 
was first suggested by L. Pitaevskii 

Further we review our results based on the stability analysis of the TDGLE for the 
rotating superfluid. The vortex nucleation in a rotating superfluid was a rather hot topic 
during 60th and 70th. Remnant vorticity did not permit to define correctly the critical 
rotational velocity for vortex nucleation at that time. We believe that next generation of 
experiments will overcome this problem, e.g., in the way suggested about 25 years ago 0. 
Together with new experimental methods of vortex detection suggested recently it can 
provide experimental test of the theory. 

Formation of vortices under a rapid quench is recognized as a fundamental problem of 
contemporary physics ||. Superfluid 3 He offers a unique "testing ground" for rapid phase 
transitions @]. Recent experiments where a rotating superfluid 3 He was locally heated well 
above the critical temperature by absorption of neutrons || revealed vortex formation under 
a rapid second-order phase transition. We will discuss the dynamics of vortex nucleation 
under a rapid quench in the framework of the TDGLE. 

II. SPIN-UP AND NUCLEATION OF VORTICES IN SUPERFLUID HELIUM 

Let us consider a cell containing a superfluid helium rather close to the superfluid tran- 
sition temperature T\. When the cell is rotated with an angular velocity Q, the normal 
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component is involved into a solid body rotation with V n = VL x r. The superfluid com- 
ponent cannot participate in the uniform rotation up to the point where the potential flow 
condition is satisfied, i.e. V x V s — 0. 

A conventional approach to the spin-up problem is to describe it by a two-fluid hydrody- 
namic model corrected by an equation for the vortex line dynamics fllOyiTl . The vortex lines 
interact with the normal component that leads to mutual friction ||11|| . It is evident that 



this hydrodynamic description does not catch the key point of the spin-up problem, namely, 
the vortex nucleation which actually causes the spin- up of the superfluid component, the 
primary superfluid relaxation mechanism toward a steady rotation. 

Another approach to the spin-up problem, which can describe both the dynamics and 
the nucleation of quantized vortices in a superfluid helium, is to use TDGLE together with 
a two-fluid hydrodynamic model [[J. Outside the vortex core which is normal, one gets the 
superfluid velocity circulation around a single quantum vortex n = 2iTh/m Correspond- 



ing set of equations, which describe the dynamics of the complex order parameter of the 
superfluid condensate, \I/ = \^/ \ exp(ix), in a rigid steady rotation of the normal component, 
looks in the scaled variables as follows @( see for details [0): 

d t ty = ~ + * - |^| 2 ^) + - ((V - iV n ) 2 ^ + ^ - |^| 2 ^) (1) 

where V n = Qr, V s = Vx, and A is a temperature-dependent parameter. 

Equation (P is reminiscent of that of the Ginzburg-Landau equation for superconductors 
in the London limit [fL3| . The role of an external magnetic field is played by the angular 



velocity, and of the corresponding vector-potential by the velocity of the normal component. 
Then by analogy one expects that at Q < Q c i, there exists a motionless superfluid component 
with no vortices. At Q > Q c i vortices will be nucleated and penetrate into the fluid producing 
a vortex lattice in the interior of a helium container. As follows from the experiments on 
the superfluid 4 He, Q c i is too low to be detected |TTJ. On the other hand, Q C 2 (which is 
analogous to H c2 in superconductors and at which superfluidity will be completely destroyed 
in the sample) is too high to be reached experimentally. 



Equation (Jl]) describes the spin-up of the superfluid part in a rigidly rotating flow of the 
normal component. 

This equation is asymptotically correct in the vicinity of the A point. Moreover, one 
can speculate that the equation qualitatively describes some aspects of spin-up for 4 He near 
zero temperature with A — > 0. Therefore, we will consider Eq. ([I]) for all temperatures with 
A -> for T -> and A ~ (T A - T)" 1 / 3 for T -> T A . 

The critical rotational velocity, Q c , for the onset of the vortex nucleation for T — > T\ 
can be found from linear stability analysis of a stationary solution. Then for this solution 
X = 0, and F — \^f\ is defined by the following equation: 

ftF+^lL + F-F*- n 2 r 2 F = 0. (2) 
r 

Equation (g) has to be complimented by the conditions at the r = and the condition at 
the outer wall r = R, where R is the radius of the container. As a boundary condition at 
the wall we take a condition of finite suppression of the superfluid density by the wall, i.e., 
d r ty + 7^ = for r = R where 7 characterizes the suppression of the order parameter. For 
7 — > we have no- flux boundary condition (d r ^> = 0). Solution of Eq. (Q) for arbitrary Q 
and R is accessible only numerically. Selected results are presented in Fig.l of Ref . 



The stationary solution is stable for Q < Q c and looses its stability above the critical 
angular velocity Q c . Instability of the stationary solution leads to nucleation of vortices 
and the corresponding spin-up of superfluid. By substitution of a solution of the form 
if = F(r) + W(r, 6, t), where W is a small generic perturbation, one obtains a linear equation 
for W. Q c is found from the existence condition of the first nontrivial eigenmode satisfying 



the boundary conditions |12j . In the limit of R 3> 1 the solution is found by matching of the 
bulk solution with the solution near the wall. The analysis shows that the most unstable 
eigenmodes are localized in the narrow layer of the width near the container wall. We 
obtained ~ \f~R <^ R for large R. The value of most unstable azimuthal number n and 
the critical frequency Q c for the container radius R is given by the expressions: 



n = Q( 7 )^ 4 , V c = yi + ^. (3) 



The dimensionless parameters Q(j), A (7) are the functions of the suppression rate 7, are 



obtained by the matching of outer and inner expansions, and are shown in Fig. 2 of Ref |T2 
One can make estimates of fl c based on Eq.(3). Indeed, in dimensional variables one gets 
Vl c = (h/mR d )[e 2 / 3 /(V3£ )] = 3.3 x 10 3 e 2 / 3 /R d 1 /sec, where e = {T\—T)/T\, £0 = 2.74x 10~ 8 
cm is the correlation length far from T\, and the power 2/3 is introduced to assure correct 
scaling of the superffuid density near the A point [^4| . The conventional Feynman equation 



| iT| has different scaling with R and no temperature dependence. On the other hand, as we 



discussed in Ref. [ 12j , the temperature dependence of Q c found is the same as in the theory 



of thermal nucleation of quantized rings |[L1]| . However, the nucleation rate in our case was 
calculated up to prefactor (see Ref. fl2fl), while in the former theory it was obtained from 
heuristic arguments. Moreover, we consider non-uniform distribution of p s due to V n that 
is impossible task for the equilibrium theory. The latter can become significant at small 
values of e and R. Thus, at e = 10~ 6 and R d = 0.1cm scaled R is of the order 10 3 , and the 
correction to the frequency from Eq.(3) can be of the order of several percent. 

We performed numerical simulations of Eq. ([!]), which details are presented in Ref [|12 



We observed nucleation and consequent tearing off of the vortices for Q > Q c irrespective 
of A. However, the character of the nucleation and asymptotic states depends on A. For 
A — > 00 and slightly above Q c we observed nucleation of several vortices. Nucleation occurs 
at nonlinear stage of the instability when a set of single zeros (four zeros for R = 65) is torn 
off at the radius R. These zeros are the seeds for the vortex cores. The vortices propagate 
into the interior of the container and finally form a perfect vortex lattice, reminiscent of 
that of the Abrikosov lattice ]13| (see Fig. 1). Further increase of Q results in formation of 
additional vortices. 

III. STABILITY OF MULTICHARGED VORTICES 

As we pointed out above, Eq.(|l|) can describe qualitatively the spin- up of superfluid 
helium also at T — > with A — > 0. Then the TDGLE becomes the nonlinear Schodinger 
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equation (NLSE). Numerical simulations show the character of the vortex nucleation is 
drastically different for A — > 0. Typically whole clusters of vortices are torn off. These 
clusters can be considered as a perturbed multicharged vortex. The multicharged vortex is 
unstable and breaks down into single charged vortices. However, the lifetime happens to be 
proportional to A and diverges for A — ► 0. The multicharged vortices with the topological 
charge ±n are known to have higher energy and decay into n single-charged, or elementary 
vortices. Since in the NLSE with small dissipation the decay time of the multicharged 
vortices can be arbitrarily large, one can expect to detect them in experiments at very low 
T. Some indirect indications of the multicharged vortices can be found in several experiments 
on vortex nucleation ||11|| . 

We consider the perturbative solution \I/ = [F(r) + rj(x, y,t)] exp[m#] of the NLSE with 



dissipation 16 



d t V = (e + i) (A^ + ^- |^| 2 ^) , (4) 

Here 77 is the complex function, and e <C 1 is the phenomenological parameter which describes 
the bulk dissipation of superflow towards the condensate. We assume that the only channel 
for the bulk dissipation at T — > is the absorption of acoustic excitations of superflow by the 
normal component which presents, e.g., due to normal 3 He atoms. This assumption is based 
on consideration that the excess energy of n-charged vortex can decay at T — > only by 
acoustic radiation of the bounding energy of n single charges. In the presence of the energy 
conservation and other integrals of motion in the NLSE this transformation of the bounding 
energy into acoustic field is a very slow process, i.e. multiple vortex can be long-lived. 
Numerical analysis of the stability problem reveals that the lifetime of the multicharged 
vortex diverges as r ~ (eAi)^ 1 and is formally infinite for NLSE. However, one has slower 
(non-exponential) instability mode at e = 0. In particularly, generic perturbations grow 
linearly in time. In this sense the multicharged vortex is metastable and may exist for a 
very long time. The instability has a nonlinear nature and originates from the interaction 
between localized eigenmodes and the continuous spectrum of the vortex radiation. The 
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crucial point here is that the simultaneous existence of the localized eigenmodes with n > 2 
and extended excitations does not contradict to the energy conservation of the Hamiltonian 
NLSE, since they contribute to the energy with opposite signs. This process is similar to the 



growth of waves with negative energy []I7[|. In Fig. 2 one clearly sees that a rotating double- 
charged and triple-charged vortices radiate away the acoustic waves. Thus, the decaying 
multicharged vortices are effective source of the acoustic radiation. 



IV. NUCLEATION OF VORTICES BY RAPID THERMAL QUENCH 

Nucleation of vortices by neutron irradiation in 3 He-B in the presence of rotation was 
studied experimentally in Ref. [§]. Ignoring non-relevant complexity of the 3 He-B specific 
multicomponent order parameter, we will use the TDGLE for a scalar order parameter ip 

w 

d t i> = Ai/; + (1 - /(r, t))iP - + C(r, t). (5) 

Close to T\ the local temperature is controlled by normal-state heat diffusion and evolves 
as f(r,t) = E exp(—r 2 /at)t~ 3 ^ 2 , where a is the normalized diffusion coefficient. E 3> 1 
determines the initial temperature of the hot bubble T* due to the nuclear reaction between 
the neutron and 3 He atom and is proportional to the deposited energy £ . The Langevin 
force ( with the correlator (CO = 2Ty5(r — r')5(t — t') describes thermal fluctuations with 
a strength Tf. 

Numerical simulations of Eq. (|5]) in 2D and 3D shows that without fluctuations (Tf = 0) 
the vortex rings nucleate upon the passage of the thermal front, Fig. 3a,b. Not all of the rings 
survive: the small ones collapse and only the big ones grow. Although the vortex lines are 
centered around the point of the quench, they exhibit a certain degree of entanglement. After 
a long transient period, most of the vortex rings reconnect and form the almost axisymmetric 
configuration. We find that the fluctuations have a strong effect at early stages: the vortices 
nucleate not only at the normal-superfluid interface, but also in the bulk of the supercooled 
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region (Fig. 3c), according to the Kibble-Zurek "cosmological" mechanism ||. However, 
later on, small vortex rings in the interior collapse, and only larger rings (primary vortices) 
survive and expand (Fig. 3d). 

We conclude that the primary source of vortices in our numerical simulations is the 
instability of normal-superffuid interface in the presence of flow. The analysis results in the 
following expression for the number of survived vortices 



where (3 = const, while v s and v c are the imposed and critical GL superflow velocity, 
respectively. This estimate is in agreement with simulations, see Ref . |§ . Eq. (||) exhibits a 
slow logarithmic dependence of the number of vortices on the level of fluctuations. For the 
experimental values of the parameters our analysis results in about 10 surviving vortices per 
heating event. It is consistent with Ref. || where as many as 6-20 vortices per neutron were 
detected. 
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FIGURES 

FIG. 1. The sequence of grey-coded images of demonstrating nucleation of vortices and 
creation of the vortex lattice. The dark shade corresponds to zero of |^|; the white one to its 
maximum value. Vortices are seen as black dots. The parameters are: f2 = 0.01, R = 65, A » 1. 
The initial condition is VP = 1 plus small amplitude broad-band noise, a) t = 100A, b) t = 200A, 
c) t = 1100A. 

FIG. 2. Number of vortices as a function of fi. a) R = 65; b) R = 140. 

FIG. 3. Gray-coded snapshots of \^>(x, y)| 2 (zero is shown in black, |^| = 1 in white) for dou- 
ble-charged vortex(a)-(c); and triple-charged (d)-(f) at the moments of time : (a) t=1700; (b) 
t=2000; (c) t=2500;(d) t=1700; (e) t=2000; and (f) t=2500. The parameters of the simulations 
are: The domain size 100 x 100 units, number of FFT harmonics 128 x 128; e = 0.001, boundary 
conditions no-flux; initial conditions were slightly perturbd double-charged vortex. Single vortices 
are presented by black spots, the acoustic field is seen in gray shade. 
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